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ABSTRACT 



We estimate the number of old, radiative supernova remnants (SNRs) detectable in 
Hi 21-cm emission line in the Galaxy. We assume that old SNRs consist of expanding 
Hi shells and that they are visible if the line-of-sight velocities are sufficiently outside 
the velocity range of the Galactic background Hi emission. This criterion of visibility 
makes it possible to calculate the background contamination and to make a comparison 
with observation. The Galactic disk in our model is filled with atomic gas of moderate 
(~ 0.1 cm~ 3 ) density representing the warm neutral interstellar medium. We assume 
that only Type la supernovae produce isolated SNRs with expanding Hi shells, or "Hi 
SNRs". According to our result, the contamination due to the Galactic background Hi 
emission limits the number of visible SNRs to ~ 270, or ~ 9% of the total Hi SNRs. 
They are concentrated along the loci of tangential points. The telescope sensitivity 
further limits the number. We compare the result with observations to find that the 
observed number 25) of Hi SNRs is much less than the expected. A plausible 
explanation is that previous observational studies, which were made toward the SNRs 
identified mostly in radio continuum, missed most of the Hi SNRs because they are 
too faint to be visible in radio continuum. We propose that the faint, extended Hi 
21-cm emission line wings protruding from the Galactic background Hi emission in 
large-scale (£, v) diagrams could be possible candidates for Hi SNRs, although our 
preliminary result shows that their number is considerably less than the expected in 
the inner Galaxy. We conclude that a possible explanation for the small number of Hi 
SNRs in the inner Galaxy is that the interstellar space there is largely filled with a 
very tenuous gas as in the three-phase ISM model, not with the warm neutral medium 
of moderate density. 
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1 INTRODUCTION 

A general picture of old SNRs is a hot bubble surrounded by 
a dense, atomic shell. The shell is preceded by SNR shock, 
which is radiative. A shock is referred to as "radiative" if the 
cooling time of the shock-heated material is shorter than the 
age of the shock. For the SNR shock propagating in a uni- 
form, homogeneous medium of hydrogen density ~ 1 cm -3 , 
it becomes radiative when the SNR is ~ 3.6 x 10 4 yrs old, or 
when the shock velocity drops to ~ 180 km s _1 (see Section 
2.3). The shell decelerates as it sweeps out more material, 
and, eventually, when the SNR is ~ 10 6 yrs old, the ve- 
locity of the shell drops to an order of the ambient sound 
speed (~ 10 km s^ 1 ) and the SNR merges into the general 
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interstellar medium (ISM). We may call old SNRs with fast- 
expanding Hi shells "Hi SNRs" because they are observable 
in Hi 21-cm emission line. 

A systematic Hi 21-cm line observation of SNRs to de- 
tect the expanding shells have been made by Koo & Heiles 
(1991; hereafter KH91). It is difficult to detect expand- 
ing Hi shells, because most known SNRs are located in 
the Galactic plane where the Galactic background Hi emis- 
sion causes severe contamination. There had been papers 
reporting the detection of Hi shells with expansion veloc- 
ities smaller than ^20 km s -1 , but the association was 
rather ambiguous (lAs^ousa 1- fc - FMces] Il973t iKnapp fc Kerrl 
Il974t IColomb fc Dubnerlll980t I1982T) . If the expansion ve- 
locity of the shell is very large, larger than the maximum 
velocity permitted by the Galactic rotation, however, then 
it could be easily discernable from the background emission. 
Fig. 1 shows such an example, where we can clearly see 
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an extended high-velocity (HV) excess emission localized at 
the position of the SNR HB 21. The positional coincidence 
and the forbidden velocity strongly suggest that the excess 
emission is emitted from the gas accelerated by the SNR 
shock. KH91 carried out a sensitive survey of Hi 21 cm emis- 
sion lines toward 103 northern Galactic SNRs and detected 
such HV gas toward 15 SNRs 1 , including 3-4 SNRs known 
prior to the survey (DeNoyer 1978; Giovanelli & Haynes 
1979; Landecker, Roger, & Higgs 1980; Braun & Strom 
1986; Koo et al. 1990). Koo, Kang, & McClure-Griffiths 
(2003) searched for similar Hi features toward 97 south- 
ern SNRs using the Southern Galactic Plane Survey data 
(McClure-Griffiths 2001) and identified another 10 SNRs. 
These 25(= 15 + 10) SNRs are candidates for Hi SNRs be- 
cause the surveys were made using single-dish radio tele- 
scopes with a large (36' or 16') beam size and high-resolution 
observational studies are necessary for confirming the nature 
of the HV emission. Such high-resolution studies have been 
done for some of the SNRs and confirmed t he shell-like struc - 
ture of the HV Hi feature s, e.g., C TB 80 jKoo et aUll99Cft . 
W44 llKoo fc Heileslll995h. W51C (iKoo fc Moonlll997h. an d 
IC 443 jGiovanelli fc Havneslll979l: iBraun fc Stromlll98rj) . 
But molecular lines studies showed that the latter three 
SNRs exhibit evidences for the interaction with molecular 
clouds, e.g., broad molecular lines, high-ratio of high- to 
low-transition CO lines, 1720-MHz maser line of the OH 
molecule, etc. (see Koo 2003 and references therein). The 
HV Hi gas in these SNRs is possibly formed by a radia- 
tive shock propagating through molecular clouds, not by the 
shock propagating through the general diffuse ISM. In sum- 
mary, including the ones interacting with molecular clouds, 
there are 25 Hi SNR candidates which appear to have shock- 
accelerated Hi gas at velocities forbidden by the Galactic 
rotation. 

How does this compare to the expected number of Hi 
SNRs in the Galaxy? If we naively multiply the total su- 
pernova rate (~ 0.02 yr" 1 ) in the Galaxy to the lifetime 
(~ 1 x 10 6 yr) of SNRs, we get ~ 2 x 10 4 . The majority of 
supernovae (SNe), however, are core-collapse SNe with mas- 
sive progenitors. Massive stars are born in associations, and 
most of these SNe dissipate their energy in producing super- 
bubbles instead of isolated SNRs. Also, a significant fraction 
of SNe occur in the Galactic halo where they merge into the 
ambient ISM before cooling becomes important. And, if the 
structure of the ISM is clo se to the three-phase model of the 
iMcKee fc Ostrikeil il977f) so that most of the interstellar 
space is filled with hot tenuous (3.5 x 10~ 3 cm -3 ) medium, 
the SNRs may overlap with other SNRs before dense shell 
formation. Meanwhile, the majority of Hi SNRs may not 
be visible because of the contamination due to the Galactic 
background emission and the limited sensitivity. The effect 
of the background contamination is difficult to calculate in 
general, which hampers the estimation of the number of vis- 
ible SNRs. But, for the Hi SNRs, if we limit to the ones 

1 KH91 observed each SNR at 9 points in a cross pattern centered 
at its catalog position. They searched for excess emission wider 
than 10 km s -1 and divided the SNRs into 3 ranks, in which 
increasing number implies increasing reliability of the detected 
Hi feature. These 15 SNRs include the fourteen rank 3 SNRs 
(excluding G117. 4+5.0 which is not considered as a SNR any 
more [Green 2001]) and the classical source IC 443. 



with large expansion velocities as in Fig. 1, it is straight- 
forward to calculate the effect because the boundary of the 
background emission in Fig. 1 is mainly determined by the 
Galactic rotation. Therefore, the statistics of such Hi SNRs 
can be related to the Galactic distribution of SNe, their en- 
vironments, and the evolution of their remnants. 

In this paper, we develop a model to estimate the vis- 
ibility of Hi SNRs in the Galaxy. The outline of the paper 
is as follows: in Section 2 we describe our model and for- 
mulation. We model the Galactic disk as an axi-symmetric, 
uniform, homogeneous disk with a central hole. We estimate 
the frequency and distribution of SNe in the disk using re- 
cent observational results on extragalactic SNe and old disk 
stars. We summarize the dynamical evolution of old SNRs 
in homogeneous, uniform medium and devise the detection 
probability. Section 3 presents the results, where we show 
the expected distribution of Hi SNRs in the Galactic plane 
and their statistics. We explore how the background contam- 
ination and the telescope sensitivity limit the visible number 
of Hi SNRs. In Section 4, we discuss the effect of embedded 
clouds and compare the result with observations. 



2 DESCRIPTION OF MODEL 

2.1 Galactic Hi Disk and Hi SNRs 

What would be the nature of the ISM that a typical iso- 
lated SNR experiences? This is perhaps the very essential 
question in deriving the number of Hi SNRs because the 
formation and evolution of dense radiative shells depend on 
the physical properties of the ambient ISM. A general pic- 
ture of the ISM is dense clouds immersed in a diffuse in- 
tercloud medium, perhaps except regions surrounding OB 
stars where the clouds might have vanished through pho- 
toevaporation and rocket effect due to strong UV radiation 
(e.g., McKee, Van Buren, & Lazareff 1984). The birthplaces 
for the isolated SNe are presumably far from OB stars, so 
as to be close to the general ISM. Then, although the em- 
bedded clouds may modify the internal structure through 
evaporation, it is the intercloud medium that primarily con- 
trols the dynamical evolution of SNRs (e.g., Cowie, McKee, 
& Ostriker 1981), and the question is what the physical con- 
ditions of the intercloud medium are. There have been sub- 
stantial studies on this subject, and it is generally consid- 
ered that the major component of the intercloud medium 
is either the warm neutral medium (WNM) with a typical 
density of ~ 0.1 cm -3 as in the two-phase model of Field, 
Goldsmith, & Habing (1969) or the hot ionized medium 
(HIM) of 3.5 x 10 ~ 3 cm" 3 as in the three-phase model of 
IMcKee fc Ostrikeil dl977ft (see McKee 1995 for a review). 
Relative filling factors of the two phases, however, are poorly 
known even in the solar neighborhood. If it is the HIM, then 
most of the SNRs might dissipate their energy without form- 
ing a radiative shell, and there will be only a few Hi SNRs. 
In the following, we calculate the number of Hi SNRs as- 
suming that the intercloud medium is composed with the 
WNM only, and, by comparing the result with observations, 
we will infer the structure of the ISM. 

We assume that the Hi gas in the Galaxy is confined to 
an axi-symmetric, center-emptied disk with an inner edge at 
3.5 kpc and outer edge at 15 kpc. The total Hi column den- 
sity perpendicular to the Galactic plane is approximately 
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constant between 3.5 kpc and 20 kpc, but has a central 
hole. The inner edge represen ts the boundary where th e Hi 
surface density drops rapidly (|Kulkarni fc Heiieslll987f l. We 
take the outer edge of the disk at 15 kpc (instead of 20 kpc) 
considering that the distribution of isolated SNe are possi- 
bly truncated at £j 15 kpc (see Section 2.2). Even if there is 
no truncation, the small disk size should not notably affect 
the result because, as we discuss in Section 2.2, the isolated 
SNe are considered to be centrally concentrated with a short 
(~ 3 kpc) exponential radial scale length. We assume that 
the thickness of the disk is constant and the density distribu- 
tion is uniform. The real Galactic disk is warped and flares 
outside the solar circle (Kulkarni, Blitz, & Heiles 1982; Lock- 
man 1984). But this complexity should have little effect on 
the result because, firstly, most (~ 80%) SNe occur within 
the solar circle and, secondly, the st ellar disk and, therefore 
the SNe disk, warps and flares too IlLopez-Corredoira et alJ 
2002). Our formulation is essentially two-dimensional and 
the primary purpose of considering the thickness of the 
disk is to estimate the fraction of SNe that explodes in the 
very tenuous halo, which may not evolve to Hi SNRs (Sec- 
tion 2.2). According to KH91, all the Hi SNR candidates 
in the inner Galaxy is within z < 200 pc, while, in the 
outer Galaxy, the one (CTA 1) located at the farthest from 
the midplane is at (r, z) = (10 kpc, 370 pc) where r is the 
Galactocentric radius and z is the height above the mid- 
plane. The latter corresponds to 314 pc at r — r© = 8.5 kpc 
if we consider that the thickness of the disk flares almost 
linearly outside the solar circle. Considering this, we adopt 
the half-thickness of our disk to be 320 pc, i.e., SNe oc- 
curring outside this disk are assumed not to form an Hi 
shell. This is greater than the half-width at half-maximum 
(265 pc) of the thick Gaussian component and less than the 
scale hei ght (403 pc) of the expo nential component of the 
Hi disk iDickev fc Lockmanl ll990'l. The total intercloud Hi 
column density perpendicular to the Galact ic pla ne in the 
solar neighborhood is ~ 90 cm -3 pc and we 

require our model disk to have the same vertical column 
density. Therefore, the density of our disk no = 0.14 cm" . 

For such axi-symmetric, two-dimensional disk, the num- 
ber of observable Hi SNRs can be written as 

N = (E SN /Hi)r/ obs , (1) 

where Esn is the total isolated SN rate in the Galaxy, /hi is 
the fraction of isolated SNe that produces Hi SNRs (i.e., the 
fraction within the center-emptied model gaseous disk) , r is 
the lifetime of Hi SNRs, and / b s the fraction of Hi SNRs 
that are detectable: 

/obs = ; / o-snM / p(x,y,i)dtdA (2) 

^sn/hit J J 

where usn(0 is the surface density distribution of SNe, i.e., 
the SN rate per unit area as a function of Galactocentric ra- 
dius r, and p(x, y, i) is the probability that a Hi SNR of age t 
is detectable considering the background contamination and 
the telescope sensitivity, (x, y) represents a coordinate grid 
centered at the Galactic center and the integration range 
over (x, y) is from r = 3.5 kpc to 15 kpc. If the detection 
probability p(x,y,t) is independent of {x,y) and is equal to 
1 during t, then / b s = 1 and TV = Esn/hit. In the follow- 
ing sections, we estimate the individual items in equations 
(1) and (2). 



2.2 Frequency and Distribution of SNe 

For our purpose, SNe are divided into two types: Supernovae 
of type la (SNe la) which have old disk stars as their pro- 
genitors and the "core-collapse" SNe (e.g., type Ib/c and 
type II SNe; hereafter SNe II) which have massive progen- 
itors. SNe II may be further divided into two groups; ones 
exploding in associations and the other exploding isolated. 
Most of the former SNe are expected to explode within a 
superbubble where the density is very low and pressure is 
very high. Such SNRs might dissipate their energy before 
radiative cooling becomes important, and we consider that 
they do not evolve to Hi SNRs. The latter group is con- 
sisted of "runaway" OB stars which are ejected from their 
paren t association and also possibly the OB stars formed iso- 
lated. [Gies^Boitonl (|l98(| derived that 10-25% of O stars 
are runaways while only about 2% of early B stars , which 
contri bute most of the SNe, are runaways. This led lMcKed 
( 1995) to conclude that there is no need to assume a sep- 
arate populati on of isolated OB stars. On the other hand, 
iFerrierel fo)95|) an alyzed the catalog of OB stars compiled 
bv lGied Jl987h and lHumphrevs fc McElrovl dl984h . and con- 
cluded that 40% of SNe II occurs isolated (see also McKee 
1990). Among the historical SNe II, at least Crab is not as- 
sociated with a grou p of OB stars and its progenitor coul d 
have been a runaway iMdzinarishvilli fc Dzigvashvillil200ll) . 
We assume that all SNe II explode in groups and that only 
SNe la produce isolated SNRs that can evolve to Hi SNRs. 
The basic idea is to choose parameters to minimize the num- 
ber of Hi SNRs when they are uncertain, so that the derived 
number is to be a lower limit. If the population of isolated 
OB stars is indeed significant, then the expected number of 
Hi SNRs will be larger than the calculated below. 

The Galactic SN la rate is inferred from the extra- 
galactic SN rates, which are usually expressed in SNu or 
SN (100 yr)" 1 (10 10 L B / X©)" 1 where L B is the blue lu- 
minosity of the galaxy (e.g., van den Bergh & McClure 
1994; Cappellaro, Evans, & Turatto 1999). Our Galaxy is 
considered to be a spiral of Hu bble type Sb or Sbc with 
L B = (2.3 ± 0.6) x 10 10 L© llvan den Berghlll98gft . and we 
take the arithmeti c mean of the average ra tes of SOa-Sb and 
Sbc-Sd galaxies in lCappellaro et alJ (Il999l) to obtain the to- 
tal SN la rate of 0.45 x 10~ 2 yr -1 . (For comparison, the total 
SN II rate is 1.8 x 10" 2 yr" 1 .) It is noted that this SN la rate 
is hi gher than the one (0.3-0.4 x 10~ 2 ) d erived from the nova 
rate llvan den Bergh fc Tammarm ^^9jJ) , but consistent with 
the es timate (0.3-0.6 xl0~ 2 ) of Ivan den Bergh fc McClurd 

SNe are centrally concentrated and vertically strati- 
fied, presumably following the distribution of their progen- 
itors. The progenitors of SN la are mostly old disk stars. 
Bulge stars could contribute as much as 15% to the to- 
tal SN la rate if the SN rate scales with the mass ratio 
of the bulge (~ 10 10 M© ) to the disk (~ 7 x 10 10 M ) 
llDawson. fc Johnsonlll994h . But the rate of SN la within 1 
kpc of the center of spiral galaxies has been found to be sig- 
nificantly lower than outside this region (Wang, Hoflich, & 
Wheeler 1997; Hatano, Branch, & Deaton 1998) , which sug- 
gests that the SN la occurs in the bulge not as frequently 
as expected. The space density of old disk stars perpen- 
dicular to the Galactic plane can be accurately represented 
by the sum of two exponentials, one with a characteristic 
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scale height of ~ 300 pc and the other which h as a scale 
height of ~ 1 kpc (e.g.. lBinnev fe MerrifieTdlll998l) . The lat- 
ter "thick disk" stars appear older ( > 10 10 yr) and more 
metal poor ([Fe/Hl ^ — .4) than the former "thin" disk 
jBinnev fc Merrifieldll998h . The total mass of the thick disk 
stars is < 10% of the thin disk stars (Trimble 2000). There- 
fore, if the SN la rates of thin and thick disks are propor- 
tional to their masses, the thick disk could contribute at 
most ^ 10% of the thin disk to the total SN la rate. We ne- 
glect the small uncertain contributions from the bulge and 
thick disk, and assume that the old stellar disk representing 
the SN la progenitors has a vertical exponential distribution 
with a scale height of 300 pc at r = r©. This implies that, 
in the solar neighborhood, 34% of SNe la explodes outside 
our ±320 pc-thick gaseous disk. We further assume that this 
fraction remains constant over the full disk, e.g., we assume 
that the fraction of SNe la exploding in the 'Hi halo' and not 
evolving to Hi SNRs /halo = 0.34. This is not unreasonable 
because t he old stellar disk warps and flares as the gaseou s 
disk does iFreudenreictJl998llL6pez-Corredoira et all2 002). 
But the fla ring of the stellar disk may st art well inside the 
solar circle llLopez-Corredoira et alj|2002fl . in which case the 
assumption of the constant fraction underestimates the SN 
la rate in the inner Galaxy. 

The radial surface distribution of old stellar disk is usu- 
ally described by a single exponential profile, but the radial 
scale lengths quoted in published papers span an enormous 
range: 1.8 to 6 kpc (Kent, Dame, & Fazio 1991 and refer- 
ences therein) . Recent near-infrared studies based on the the 
COBEfDIRBE data or the Two Micron All Sky Survey data 
appear to yield a shorter scale length (2.1-2.8 kpc) consis- 
tently (Freudenreich 1998; Ojha 2001; Drimmel & Spergel 
2001 and references therein). On the other hand, studies on 
extragalactic SNe show that, although the number of SN la 
is small (54) to constrain the slope very well, the distribu- 
tion could be represented by a single radial scale length of 
^25/5.5 where R25 is the semi major axes of isophot e hav- 
ing 25.0 mag arcsec -2 in blue llvan den Bergbl ll997). The 
Andromeda galaxy (M31), for example, has R25 = 20 kpc 
(Trimble 2000) , so that the radial scale length of SN la would 
be 3.6 kpc. We adopt 3.0 kpc as the radial scale length of the 
SN la surface density distribution in the Galaxy, and assume 
that the distribution extends all the way to the Galactic cen- 
ter. The stellar disk may have a ce ntral hole where a stellar 
bar is located iFreudenreicbl ll998) . But, since the surface 
density profile averaged over the azimuth remains ap prox- 
imately exponential to the center iFreudenreicbllT998T) . the 
assumption of an axisymmetric exponential distribution is 
plausible. Hence, after normalization, the SN rate per unit 
area in our model disk is given by 



csnM = 00 exp(— r/3.0 kpc) 



(3) 



where ao = 5.5 x 10 -5 kpc -2 yr -1 . With the above distri- 
bution, the fraction exploding inside the central hole /hole = 
0.34. Therefore, the fraction of isolated SNe that explodes 
within our model gaseous disk /hi = (1 — /halo)(l — /hole) = 
0.44. The frequency of Hi SNRs in our model Hi disk, there- 
fore, is Ssn/hi = 0.20 x 10~ 2 yr -1 , or one per 500 yrs. 



2.3 Dynamical Evolution of Hi SNRs 

In our model, SNRs evolve in a homogeneous, uniform 
medium and we are interested in their late-stage evolution. 
(The effect of embedded clouds will be discussed in Sec- 
tion 4.1.) SNRs form a dense, neutral shell through radia- 
tive cooling when the cooling time of the shock-heated ma- 
terial becomes shorter than the age of the SNRs. Numer- 
ical studies showed that the shell formation occurs 'catas- 
trophically', accompa nying multiple shocks and oscillations 
of shock velocity (e .g., lFallelll98ll iKimoto fc Chernofj|l997l : 
iBlondin et al.lll998l) . The shell is also dynamically un stable 
to develop nonspherical structure felondin et alj[l 998). But 
this violent transition lasts only for a relatively short time, 
and the global dynamics of radiative SNRs for most of time 
could be described by a spherically-symmetric, stea dy model 
fe.g.. lCoxlll972l;IChevalier|ll974l ; lciofE et al]ll98lgh . We use 
the results of TcioilT^tHnT i r988F in the following to describe 
the SNR evolution. Their model neglects magnetic field and 
thermal conduction. Magnetic field inhibits compression to 
make the shell thicker and the outer radius larger. Thermal 
conduction lowers the temperature and increases the den- 
sity in the hot bubble, the effectiveness of which depends 
on the magnetic field configuration. These effects, however, 
are dynamically important only at very late times when the 
shell is sufficiently slowed down. For example, the radius and 
veloci ty of the shell in the numerical result of ISlavin fc Coxl 
(1992), which includes a uniform 5 fiG magnetic field and 
unimpeded thermal conduction, are described well by the 
formulae in this section until the expansion velocity drops to 
~ 50 km s _1 . For a tangled random magnetic field, the shell 
will be more compres sible and the effect of the conduction 
will be reduced (e.g.. lMcKeelll995l) . so that the agreement 
will persist at lower velocities. Since the visible Hi SNRs 
have large (^ 50 km s" 1 ) expansion velocities (Section 2.4), 
the magnetic field and thermal conduction may be neglected 
in our formulation. 

We consider the time when the first element of shocked- 
heated gas cools to zero as the onset of the formation of Hi 
shell. This shell-formation time t s { is given by 



t Bt = 3.61 x 10 4 n - 4/7 B 5 3 1 /14 



yr, 



(4) 



where S51 is the SN energy released to the ISM in units 
of 10 5 ergs, and we assumed that the metallicity is given 
by solar abundances. After some finite violent transition pe- 
riod, the shell ex pands steadily, driving a radiative shock. 
ICioffi et alJ <ll988T) introduced an "offset" power law that can 
accurately describe the radius and velocity of the shock dur- 
ing 0.4t a f ^t^j 13t s {. They introduced a parameter (£ p d s = 
t B f/e where e is the base of natural logarithm) representing 
the beginning of the radiative, or pressure-driven snowplow 
stage when the dynamics begins to deviate from the Sedov- 
Taylor solution due to radiative losses, and normalized all 
the variables using this parameter. Instead of introducing 
this extra parameter, we slightly modify their offset power- 
law (equation 3.32) using t s { defined in equation (4) to ob- 
tain 



/ll J 1_\ 3 / 10 

UOtrf 10 J 

(11 t 1 

Vsl \wu ~ To 



(5) 
(6) 



Visibility of Old Supernova Remnants in Hi 21-cm Emission Line 5 



where the radius and velocity of the SNR shock at the shell- 
formation time are 



R st = 20.0n" 3/7 EH 7 pc, 
v sf = 179nJ /7 £~ 1/14 kms -1 . 



(7) 
(8) 



We use equations (5) and (6) to describe the evolution of the 
radius and velocity of the Hi shell from t s { to ~ 13i s f- The 
maximum disagreement between the expressions in equation 
(5) and the original expression of Cioffi et al. is < 0.1% 
during this time period. For E51 = 1 and no = 0.14 cm -3 , 
we have t s f = 1.11 x 10 5 yr, R s { = 46.4 pc, and v a f = 135 km 
s- 1 . 

The maximum age of Hi SNRs is defined as the time 
when the shock velocity drops to about the ambient sound 
speed (~ 10 km s" 1 ) so that the interior pressure is com- 
parable t o the ambient pres sure and the shell breaks up. 
Following Cio ffi et alJ J198 8D. we define this merging time 
imerge as the time when the shock velocity drops to (3(= 2) 
times the ambient isothermal sounds speed c s . Then, us- 
ing equation (6) and writing the fractions as decimals, the 
lifetime of Hi SNRs r 



Emerge — t s f becomes 



10 

11 



Vsf 

(3c a 



10/7 



0) 



2.02 x lO 6 ^)- 1 - 43 ^- - 367 ^ 316 yr-0.91* sf) 



where c a ,a = c s /(10 6 cms -1 ). If v s t/f3c a 2> 1, the second 
term in the second expression may be dropped. Since the 
offset power-law for v a (equation 6) is accurate until t 13t s f 
or v s H 0.16« s f, equation (9) becomes less accurate for small 
( < 2) (3. For (3 = 2, c s , 6 = 1, and n = 0.14 cm" 3 , r ~ 
1.44 x 10 6 yr. Note that, since we require the expansion 
velocity of the SNR to be considerably greater than f3c a in 
order to avoid the background confusion, the observerable 
period of Hi SNRs is considerably shorter than r (Section 
2.4). 

In order to compute the detection probability in Section 
2.4, we need the mass of the shell as a function of its velocity. 
At the time of shell-formation, the swept-up mass of H atoms 
using _R s f in equation (7) is 



M Bf = 829n~ 2/7 E 6 5 { 7 M . 



(10) 



During the shell-formation stage, a good fraction ( ~ 50%) of 
this mass is compressed into a den se, neutral shell llChevalierl 
Il974l ; iMansfield fc Salpeterlll974T) . and, after the formation, 
the subsequently shocked ambient gas is simply added to 
the shell. The hot interior gas also continuously runs into 
the dense shell to add extra mass, but its contribution is 
relatively small. We assume that the shell appears with a 
mass of 5M B { at t B f and the newly swept-up mass is simply 
added to the shell. Then, from equations (5) and (6), the Hi 
mass of the shell M(v s ) is given by 



M{v a ) = M s{ 



(-) 



-9/7 



-(1-5) 



(11) 



where the second term represents the fraction of the mass 
left in hot interior at the stage of the shell formation. We 
adopt S — 0.5 following the numerical results. 



2.4 Fraction of Observable SNRs 

It is convenient to rewrite the fraction of observable Hi SNRs 
(equation 2) as an integral over velocity because the detec- 
tion probability is a function of v„: 



fo 



TSN(r) 



pjx,y,v a 
\dvs/dt\ 



■dv a dA 



where, from equation (6), 



dv s 



dt 



77 v s f 
100 t 



- ( — \ 

f \V s i) 



17/7 



(12) 



(13) 



The upper limit i> max = v B f and the lower limit u m i n = f3c s . 
If p = 1 in this velocity interval, or equivalently from i s f to 
t sf + r, then / obs = 1. 

The detection probablity p(x,y,v s ) may be expressed 
as a product of two probabilities; ps(x, y, v s ), the probabil- 
ity of detection that is limited by the confusion with the 
Galactic background Hi emission, and pT(x,y,v a ), that the 
one is limited by the telescope sensitivity. Our intention is 
to compare our results with the observations described in 
Section 1, and we define pb(x, y, v a ) so that a Hi SNR is ob- 
servable if the line-of-sight velocity of its shell is sufficiently 
outside the velocity range of the background emission, i.e., 

{1 if V SNR + V S ^ «J,max + A« cr 

Or VSNK — V s ^ U/,mln ~ At) a 
otherwise, 

where vsnr is the systematic LSR velocity of the SNR and 
ff.max and w^mm are the maximum and minimum LSR ve- 
locities toward £ permitted by the Galactic rotation. We 
compute the LSR velocities using a flat rotation curve with 
vq = 8.5 kpc and vq = 220 km s _1 . The Galactic back- 
ground Hi emission extends beyond W£ imax and u«, m i n be- 
cause of thermal motions, turbulent motions, non-circular 
motions associated with spiral density waves, etc, and Aw cr 
represents this extra velocity extent of the background emis- 
sion. Inspection of large-scale position velocity diagrams 
suggests that Av CI varies with Galactic longitudes and lat- 
itudes. We adopt Av cv = 50 km s _1 , which is sufficiently 
large in most parts of the Galaxy. The dashed lines in Fig. 1 
(and Fig. 8) show the boundary of the background emission 
defined in this way. It is obvious that the expansion veloc- 
ity of Hi SNR must be greater than 50 km s _1 in order to 
be visible. When its expansion velocity is 50 km s _1 , the 
age and radius of the Hi SNR t = 3.85i sf = 4.27 x 10 5 yr 
and R 3 = 1.537? s f = 71.0 pc from equations (5) and (6), 
respectively. 

Pt, the detection probability limited by the telescope 
sensitivity, depends on how the SNR couples with the tele- 
scope beam, i.e., whether the SNR is resolved or unresolved 
by the telescope beam. We consider two limiting cases: 

Model A: Flux-unlimited (pr = 1) — This model applies 
to an ideal case in which the visibility is not limited by the 
telescope sensitivity. The detection is only limited by the 
confusion due to the Galactic background emission. 

Model B: Flux-limtted (pi ^ 1 ) — This model applies to 
observations with limited telescope sensitivity. We consider 
only the case in which the SNRs are unresolved by the tele- 
scope beam. If a source is unresolved, the 21-cm line flux is 
proportional to the Hi mass of the source and inversely pro- 
portional to the square of the distance. Therefore, there is a 
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maximum distance ri max that it is detectable, i.e., pt = if 
d > d max . The maximum distance is derived below. 

Since the Hi 21-cm line emission from the shell is opti- 
cally thin, the integrated line intensity (expressed in antenna 
temperature Ta) is proportional to the total Hi mass M(v s ): 



T A dv = 



3hcA 10 A cB M(v B ) 
32-Kkmu d 2 



(14) 



where Aiq = 2.86888 x 1(F 15 s _1 is the spontaneous tran- 
sition probability between the hyperfine structure levels of 
Hi, A e g is the effective area of the telescope, and the other 
symbols have their usual meanings. In deriving equation 
(14), we assumed that the SNR size is much smaller than 
the telescope beam. The SNR size at the shell-formation is 
2R s{ = 92.8 pc = 32' x (d/W kpc)" 1 , and larger at later 
times. Hence, it appears that most Hi SNRs should have 
sizes larger than the beam size (~ 30') of a 25-m telescope 
(cf. Fig. 5) . However, the excess emission at extreme veloci- 
ties that we would detect comes from small cap portions of 
expanding shells, and Ta from equation (14) might yield an 
approximate average intensity of the emission. If the shell is 
thin and expands with a uniform velocity v s , then the line 
profile will be a rectangular shape with velocity extent from 
~v a to +v s , so that the left-hand side of equation (14) be- 
comes Ta x 2v„. Hence, if the minimum detectable antenna 
temperature of the observation is ATa, then from equations 
(11) and (14), the maximum distance d m ax(us) is given by 



fin 



VS \VsiJ \VstJ 



-,1/2 



(15) 



where d s { is the maximum distance at the shell-formation 
stage: 



rfsf = 



ShcAio A e fi SM st 



327rfcra H AT A 2v aC 



1/2 



(16) 



i on! A lArn \l/2rl/2 -3/14^11/28, 

= 1.89( J 4 C ff i 6/ATA,-i) ' S ' n 1 E 51 ' kpc, 

with A e fi,6 = A cff /(10 6 cm 2 ) and AT A ,-i = AT A /(0.1 K). 
The probability pt in Model B is given by 



pr(x,y,v s 



if d < dmax(«s) 

otherwise. 



For a 25-m telescope with A e s^ = 3.0 and ATa, 
dsf = 3.5 kpc for no = 0.14 cm and S = 0.5. 



= 1.0, 



3 RESULTS 

The basic results are shown in Fig. 2 where we plot the ex- 
pected surface density distribution of Hi SNRs in the Galac- 
tic plane. Fig. 2(a) shows the distribution of all Hi SNRs, 
i.e., when pb = Pt = 1. The total number of Hi SNRs in our 
model gaseous disk is Esn/hit = 2850. This is 44% of the 
total (6480) isolated SNRs in the Galaxy. The rest are in the 
halo (34% or 2200) and in the central hole (22% or 1430). 
The distribution in Fig. 2(a) is simply due to the exponential 
distribution of the SN surface density. Fig. 2(b) shows the 
distribution of observable Hi SNRs in Model A, i.e., when 
Pt = 1. The total number of observable SNRs is now 267, 
or 9.4% of the total Hi SNRs. They are concentrated along 
the loci of tangential points because the systematic veloc- 
ities of the SNRs in those regions are close to either the 



maximum or the minimum velocities of the Galactic back- 
ground emission. Only the receding (approaching) portion 
of the shell will be detectable for the Hi SNRs in the first 
(fourth) quadrant. The concentration near I = 0° is because 
the minimum (maximum) velocity permitted by the Galac- 
tic rotation is close to zero in the first (fourth) quadrant. 
For those SNRs, only the approaching (receding) portion of 
the shell will be detectable in the first (fourth) quadrant. 
Near £ = 180°, both the receding and approaching portions 
are detectable (cf. Fig. 4). Fig. 2(c) is for Model B when the 
observation is made with a radio telescope with an effective 
area of 3 x 10 6 cm 2 , which corresponds to a telescope with 
a diameter of 25 m, and ATa,-i=1.0. Now, the SNRs on 
the far side of the Galaxy are not observable and the total 
number is 96, or 3.4% of the total Hi SNRs. 

Fig. 3 (left-hand panel) is the distribution of expansion 
velocities of Hi SNRs. The distribution of all Hi SNRs is 
given by dN/dv s oc v a 11 (equation 13). Among these Hi 
SNRs, the ones with small expansion velocities are largely 
unobservable because of the contamination due to the Galac- 
tic background emission, which is shown by the Model A 
distribution. On the other hand, the Model B distribution 
shows that the SNRs with large expansion velocities are 
largely unobservable mainly because of the limited telescope 
sensitivity. Fig. 3 (right-hand pannel) is the distribution of 
distances to Hi SNRs. The bowl between d — 5 and 12 kpc 
in the distribution of all Hi SNRs is due to the central hole. 
The rapid drop beyond 7-8 kpc in Model A distribution is 
due to the small number of observable SNRs beyond tan- 
gential points in the inner Galaxy (cf. Fig. 2b). The Model 
B distribution shows that, with a 25-m telescope, we can 
detect all the observable Hi SNRs to 4 kpc and most up to 
7 kpc, but no Hi SNRs beyond 9 kpc. 

A statistical property that can be directly compared 
with the observations is the Galactic longitude distribution 
of Hi SNRs (Fig. 4). The distribution of all Hi SNRs peaks 
toward I = 0°, but not strongly because of the central hole. 
The background contamination brings down the real distri- 
bution by a factor of 5-16 depending on £ (Model A). The 
effect is largest at I = 10°-30° (or 340°-350°) and smallest 
at t — 160°-200°. The telescope sensitivity further decreases 
the number of SNRs, most notably between I — 330° and 
30° (Model B), which is understandable from Fig. 2(c). For 
Model B, Fig. 4 also shows whether the observable portions 
of the SNRs are approaching, receding, or both. Another 
statistical property that can be directly compared with the 
observations is the angular size distribution (Fig. 5). The 
distribution of all Hi SNRs peaks at 40'-50'. In the Model A 
distribution, there is a strong peak at 20-30' and a plateau 
between 30' and 70'. Because of the telescope sensitivity, 
there is no observable Hi SNRs smaller than 40' and the 
distribution peaks at 60'-70' (Model B). 



4 DISCUSSION 
4.1 Effect of Clouds 

The real ISM is pervaded by cold Hi clouds and their ef- 
fect on our results should be discussed. The evolution of 
SNRs in a cloudy ISM has been studied numerically by 
ICowie et alJ <ll98lf) . One of their models (model 6) is for the 
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two-phase ISM and applicable to our case. They adopted 
no = 0.2 cm -3 and Ebi = 0.31, and ass umed that the 
Spitzer's standard Hi clouds (jSnitzerl ll97Sl 'l are randomly 
distributed with a relatively large (7%) volume filling factor 
(see below). They simulated the evolution until the begin- 
ning of radiative phase and showed that neither cloud evap- 
oration nor the dynamical effects of the clouds affect the 
evolution of SNRs in a significant way. Their shell-formation 
time and radius agree with those from equations (4) and (7). 
They did not derive the fraction 8 of mass compressed into 
neutral shell at the time of shell formation, but, since the 
interior structure is only slightly modified, we expect that 
5 is not significantly different from the homogeneous case, 
i.e., 0.5. 

In radiative phase, the cloud-shell interaction depends 
on the ratio of their column densities. The column density of 
the shell of the visible Hi SNRs is (3.3-8.8) x 10 18 cm" 2 (cf. 
equation 11), which is less than that of typical clouds (~ lx 
10 20 cm -2 ). T hus, the collisions with clouds will punch holes 
in th e shell jMcKee fc OstrikeJ Il977t lOstriker fc McKed 
1988). The hot interior gas, however, will drive a shock 
to reform the missing portion of the shell. Since the col- 
umn density for a stro ng shock to be rad iative iV ra d ~ 
10 17 ' 5 (w s /100 kms" 1 ) 4 jMcKee et al.lll987l) . the required 
distance for the reformation of the shell is N rax i/no ~ 
0.73(« s /100 km s -1 ) 4 pc. This is much less than the mean 
distance for a particular portion of the shell to experience 
successive collisions (see below). Hence, in the two-phase 
ISM, the Hi SNRs will be punched by cold clouds as they 
expand, but the holes will be quickly recovered. The shell 
loses mass and momentum by the collision with clouds. 

The evolution of such r adiative SNRs in cloud y 
medium has been analyzed by lOstriker fc McKeel {l988). 
The key factor that determines the importance of clouds 
is exp[— (R s — i? s f)/A c i] where A c i is the cloud mean free 
path. For example, the mass in the shell at the time of 
shell formation decreases by this factor due to the collisions 
with clouds as it expands. If (R 3 — R B {)/\ C \ <C 1, there- 
fore, the 'cloud-punching' is not important. If the clouds 
are spherical with radius a c i and their volume filling fac- 
tor is / c i, then A c i = (4/3)a c i// c i. The filling factor of 
cold clouds can be inferred from their characteristic and 
space-averaged densities. The two-phase equilibrium yields 
the density of col d clouds n c i ~ 40 cm" 3 , which agrees 
with observations <lMcKedll995l) . The space-avera ged den- 
sity o f cold Hi near the Sun is ~ 0.3 cm -3 fe.g.. iFerrierd 
Il995l) . so that / d ~ 0.3/40 = 0.75%. The characteristic 
radius of the clouds inferred from the Hi c olumn density 
(N c i = 4 x 10 20 cm" 2 ) of 'standard' Hi clouds JSpitzeill978ri 
is a c i = (3/4)iVci/n c i = 2.4 pc. Therefore, A c i — 430 pc and 
(-R s --Rsf)/Aci < 5.7 x 10" 2 for visible Hi SNRs. This implies 
that the cloud punching effect is not important for the evo- 
lution of visibl e Hi SNRs and therefor e for their statistics. 

Recently iHeiles fc Trolandl J2003I) showed that in two 
regions they studied the cold ISM is distributed in huge 
blobby sheets of thickness ~ 0.11 and £ 3.6 pc with lengh- 
to-thickness aspect ratios ~ 280 and ~ 70, and, proposed 
that, if these characteristics are general, the cold ISM could 
be organized into a small number of large sheetlike struc- 
tures instead of a large number of randomly distributed 
small clumps. If the ISM morphology is close to this clumpy 
sheet model, then the probability for a SNR to cross the 



sheets will be small and only a portion of the SNR will be 
affected by the interaction, so that the the statistics of Hi 
SNRs will be determined by the nature of the intercloud 
medium. 



4.2 Comparison with Observation 

Our results can be compared with results of KH91 and 
iKoo et alJ feOOSft introduced in Section 1. In short, KH91 
did a sensitive survey of Hi 21 cm emission lines toward 103 
northern Galactic SNRs using the Hat Creek 26-m telescope 
(A e fi ~ 3.0), and detected faint extended wing at forbidden 
velocities toward 15 S NRs. 2 The sensit ivity (1.5<r) of the 
survey was ~ 0.06 K. IKoo et all (|2003l ) searched for simi- 
lar Hi features toward 97 southern SNRs and identified an- 
other 10 SNRs. The Southern Galactic Plane Survey data 
was obtained with the Parkes 64-m telescope (A c g^ — 18). 
The sensitivity (1.5cr) of the survey was 0.09-0.14 K. These 
25 (=15+10) SNRs constitute the Hi SNR candidates (see 
comment in Section 1). On the other hand, the expected 
number from Model B for the parameters of KH91, e.g., 
A eB , 6 /ATA -i = 3.0/0.6 = 5, is ~ 120, and for those of 
IKoo et al.l ( 120031) is ~190-220. Therefore, the observed num- 
ber of Hi SNRs is less than the expected by a factor of £ 5. 
Fig. 6 shows the Galactic longitude distribution of the Hi 
SNR candidates. For comparison, we show the expected dis- 
tributions and also the distribution of all known Galactic 
SNRs in the Green's catalog (Green 2001). Most candidates 
in the first quadrant have an excess emission at positive 
velocities, while those in the fourth quadrant at negative 
velocities. This seems to be consistent with our expectation 
(Fig. 4). Fig. 7 shows the angular size distribution. Most 
(72%) Galactic SNRs have sizes smaller than 30', while the 
average size of Hi SNR candidates is 46'. A significant (36%) 
fraction of Hi SNR candidates has size less than 30' whereas 
none is expected from the model. 

The discrepancy between the observed and the expected 
numbers is significant and needs to be explained. A plausible 
explanation is that the old SNRs are too faint to be iden- 
tifie d in radio contin uum, so that the Hi surveys by KH91 
and lKoo et alJ i2003l) . which were based on the current (ra- 
dio continuum-based) SNR catalog, might have missed most 
of the SNRs with Hi shells. This thought led us to search 
for Hi SNRs based on the Hi data. We used the Leiden- 
Dwingeloo (LP) Hi survey data, which was obtained by 
lHartmann fc Burtonl lll997l) using the Dwingeloo 25-m ra- 
dio telescope (HPBW=36'). The survey covers the Galac- 
tic plane (|6| Z 10°) between I = 0° and 260°, and the 
rms sensitivity of the final data cube is 0.07 K. We first 
examined about 130 SNRs in the Green's catalogue and 
have found high-velocity features toward ^ 10 SNRs. These 
SNRs look protruding from their surroundings in (£, v) dia- 
gram (cf. Fig. 1). We then examined individual (t, v) maps 
with the naked eye, and have identified ~ 70 similar-looking, 
forbidden-velocity features not related to the known SNRs. 

2 After the KH91 survey in 1991, more SNRs have been identified 
in radio continuum and in X-ray, and the number of northern 
(<5 > - 38. °5) SNRs increased from 103 to 152 (Green 2001), so 
that the number would increase to 15 x (152/103) = 22 by scaling. 
But this would not affect the discussion in this paper. 
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Fig. 8 is an example where we can see several such features. 
Some of these "forbidden-velocity (FV) wings" are associ- 
ated with galaxies or perhaps HII regions, but most (70%) 
of them are not. Could they be the so-far undiscovered old 
Hi SNRs? We have been looking into these regions in var- 
ious wavelengths, but have not identified their natures yet. 
But the distribution of FV wings appears to be very differ- 
ent from the expected distribution of visible Hi SNRs (Fig. 
9): Between I = 0° and 70°, the number of FV wings is 
much less than the expected number of Hi SNRs, while, be- 
tween i = 70° and 250° , it is significantly greater in general. 
Possible candidates for such extended forbidden-velocity Hi 
wings other than SNRs are stellar wind- blown shells, neutra l 
stellar winds associated with protostars llLizano et al Il988ft . 
high- velocity and/or intermediate clouds, galaxies, et c. The 
detailed results are in preparation jKang et alJl2003l) . 

If the observed number (^ 25) turns out to be close to 
the total number of Hi SNRs in the Galaxy, then the dis- 
crepancy implies that either some of our input parameters 
are inaccurate or some of our basic assumptions are wrong. 
The parameters and assumptions regarding the frequency 
and distribution of SNe, however, are chosen to minimize 
the number of Hi SNRs when they are uncertain, e.g., we 
have assumed that all core-collapse SNe explode in groups 
and only SNe la produce isolated SNRs. The total SN la rate 
cannot be much lower than what we adopted (0.45 x 10 -2 ) 
if we consider that there are at least two SN la (Tycho and 
SN 1006) during last 2,000 years in the solar neighborhood. 
We have discussed the effect of magnetic fields, thermal con- 
duction, and embedded clouds on the evolution of SNRs in 
Sections 2.3 and 4.1, and concluded that they are not dy- 
namically important for the visible Hi SNRs which are mov- 
ing fast. If the Hi disk is denser and thinner, the number of 
visible Hi SNRs will decrease, but not by a large factor. For 
example, if the density of the Hi disk is ~ 0.5 cm' 3 and 
its thickness is ~ 200 pc, the number decreases by ~ 40%. 
One possibility is that the interstellar space in the inner 
Galaxy is not filled with WNM, but mostly filled with a 
very hot, diffu se (3.5 x 10~ 3 cm" 3 ) med ium as in the three 
phase model of lMcKee fc OstrikeJ il977t) . In such a tenuous 
medium, the SNRs develop dense shells when they are very 
old (9 x 10 5 yr, equation 4), at which time their radius is 230 
pc. But this radius is very large, so that it is possible that the 
SNRs come into equilibrium with the interstellar pressure 
or overlap with other SNRs before den se shell formation, 
in wh ich case Hi shells are not formed jMcKee fc Ostriked 
Il977l) . Embedded clouds could be accelerated by the SNR 
shock a nd also by th e ram pressure, but numerical simula- 
tions bv lCowie et al.l dl98 J) showed that they are substan- 
tially destroyed by evaporation before acquiring high veloc- 
ities. In the three-phase ISM, the observed Hi SNRs may be 
the ones in special environments, e.g., near large atomic or 
molecular clouds. 

The above conclusion on the structure of the ISM in 
the inner Galaxy would not be totally changed even if the 
FV wing sources turn out to be Hi SNRs because, accord- 
ing to Fig. 9, there are much less number of FV wings than 
what the model predicts in the inner Galaxy. In the outer 
Galaxy, however, the number is significantly greater than 
the predicted. The model distribution can be made to have 
more Hi SNRs in the outer Galaxy, for example, by increas- 
ing the radial scale length of SNe distribution, but it cannot 



be made to be consistent with the observed FV wing dis- 
tribution simply by changing some of its parameters. We 
postpone further discussion on the FV wings and their sig- 
nificance after identifying their natures. 



5 SUMMARY 

We have developed a simple model to estimate the number 
of old, radiative SNRs visible in Hi 21-cm emission line in 
the Galaxy. Old SNRs are expected to be surrounded by 
rapidly expanding Hi shells, but the contamination due to 
the Galactic background Hi emission severely limits their 
visibility. The visible SNRs are defined as the ones with a 
sufficiently large expansion velocity, so that their line-of- 
sight velocities are significantly (^50 km s -1 ) outside the 
velocity range of the Galactic background Hi emission. An- 
other factor that determines the visibility is the telescope 
sensitivity, which is included in our formulation. Some of 
the essential features of our model are (1) the Hi gas in the 
Galaxy is confined to an axi-symmetric, center-emptied disk 
with an inner edge at 3.5 kpc and outer edge at 15 kpc, (2) 
the Hi disk is uniform and homogeneous, and the density of 
the disk is 0.14 cm -3 , (3) the disk is rotating with a constant 
speed vq = 220 km s _1 , (4) only SNe la produce isolated 
SNRs, (5) total Galactic SN la rate is 0.45 x 10~ 2 yr -1 , but 
only 44% of them explodes within the gaseous disk to pro- 
duce Hi SNRs (34% explodes above and below the disk while 
22% within the central hole), and (6) the SNe la are expo- 
nentially distributed with a radia l scale length of 3 kpc. We 
use the equations of ICioffi et al.l (^.988) in a slightly modi- 
fied form to calculate how the radius, velocity, and mass of a 
SNR shell vary with time, and, at each position in the disk, 
we derive the number of visible SNRs. 

According to our result, the number of Hi SNRs in the 
Galaxy is 2850. Galactic background contamination limits 
the number of visible SNRs to ~ 270, or ~ 9% of the total 
Hi SNRs. They are concentrated along the loci of tangential 
points and near £ — 0° (Fig. 2b). The telescope sensitiv- 
ity may prevent observing the ones on the far side of the 
Galaxy (Fig. 2c). We have calculated the expected number 
of visible Hi SNRs for realistic observational parameters and 
compared the results with observations. There have been 
two systematic Hi 21-cm emission line studies toward 103 
northern and 97 southern SNRs, and faint extended wings 
at forbi dden velocities ha ve been detected toward 25 SNRs 
(KH91; iKoo et alJl2003l) . The expected number is much 
( ^ 5) greater than this. 

We consider that our estimate is on the conservative 
side and the discrepancy is significant. A plausible explana- 
tion is that previous observational studies, which were made 
toward the SNRs identified in radio continuum, missed most 
of the Hi SNRs because they are too faint to be visible in ra- 
dio continuum. We propose that the faint extended wings at 
forbidden velocities seen in Hi 21-cm line surveys (e.g., Fig. 
8) could be possible candidates for such old Hi SNRs. Ac- 
cording to our our preliminary result, however, their Galac- 
tic longitude distribution is quite different from the expected 
distribution of visible Hi SNRs, e.g., the number of FV wings 
is considerably less than the expected in the inner Galaxy. 
The natures of the FV wings need to be identified. If the 
observed number (^ 25) turns out to be close to the total 
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number of Hi SNRs in the Galaxy, then the discrepancy im- 
plies that either some of our input parameters are inaccurate 
or some of our basic assumptions are wrong. A possible ex- 
planation is that the interstellar space in the inner Galaxy is 
largely filled with a very tenuous gas as in the three-phase 
ISM model. In this case, the observed Hi SNRs are likely 
the ones in special environments, e.g., near large atomic or 
molecular clouds. 
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Figure 1. {£, v) diagram of Hi 21-cm emission at b = 5.°0. 
Notice the excess emission at the position of the SNR HB 21 
(t = 89°) between 40 and 100 km s _1 . The dashed lines represent 
the boundaries of the Galactic background emission used in our 
model (see Section 2.4). The map has been made from th e Leiden- 
Dwingeloo Hi survey data iHartmann fc B urton 199^), and the 
contour levels start at 0.1 K and end at 10 K. 



12 B.-C. Koo and J.-h. Kang 




All 



I , -I--:- C .1. i , I , , L_ 



■J- fMM 
S. / / 

i m 


j, ii' — ^ 
'' , 




■; \s\ \ 










i L 




fy 1 


\ 

\ 






(b) V. 




Model A ' 





(c) 



RfJI 



Model B 



is la 



•> 1U IS 



15 1U 



1U 15 



Figure 2. Expected surface density (kpc -2 ) distribution of (a) all Hi SNRs (pe = PT = 1), (b) observable Hi SNRs in Model A 
(pT = 1): and (c) observable Hi SNRs in Model B (A e g ^/ ATa, — i = 3.0). The inner and outer circles marked by dashed lines are at 
r = 3.5 kpc and 8.5 kpc, respectively. In our model, the region inside the inner circle is devoid of the ISM and there the SNe explosion 
do not produce Hi SNRs. 
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Figure 3. Expected distribution of Hi SNRs in: (left-hand panel) expansion velocities; (right-hand panel) distances. 
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Figure 4. Expected Galactic longitude distribution of Hi SNRs. 
For the distribution in Model B, we mark whether the observable 
portions of the SNRs are approaching, receding, or both. 




Figure 5. Expected angular size distribution of Hi SNRs. 
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Figure 6. Galactic longitude distribution of SNRs. The distri- 
bution of the observed Hi SNR candidates is shown by a thick 
histogram, where we mark whether their approaching or receding 
portions are detected. For comparison, the expected distributions 
for the parameters of KH91 and Koo et al. (2003) are shown by 
solid and dotted lines, respectively. The distribution of all known 
231 SNRs is also shown as a thin histogram. 
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Figure 7. Angular size distribution of SNRs. The meanings 
of curves and histograms are same as those of Fig. 6. For the 
distribution of the Hi SNR candidates, wc mark the contributions 
from the northern (Hat Creek 26m) and southern (Parkcs 64m) 
sky surveys, separately. 
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Figure 8. (£, v) diagram of Hi 21-cm emission at b = — 0.°5. 
Contour levels start at 0.1 K and end at 10 K. Note the spatially- 
confined, faint extended wings at forbidden velocities (FV wings). 
Some of them are associated with SNRs, e.g., the one between 
v LSR = 80 and 140 km s" 1 near I = 49° (SNR W51C) and also, 
possibly, the one between i>lsr = 80 and 110 km s~ 1 near i = 54° 
(SNR HC 40). But the others arc not related to the known SNRs. 
The dashed lines are the same as those in Fig. 1. 
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Figure 9. Galactic longitude distribution of FV wings compared 
to the expected distribution of Hi SNRs. 



